Non-analytic dependence of the transition temperature of the homogeneous dilute 
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We show that the shift in the transition temperature of the dilute homogeneous Bose gas is 
non-analytic in the scattering amplitude, a. The first correction beyond the positive linear shift 
in a is negative and of order In a. This non-universal non-analytic structure indicates how the 
discrepancies between numerical calculations at finite a can be reconciled with calculations of the 
limit a 0, since the linearity is apparent only for anomalously small a. 
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The shift of the critical temperature of the dilute ho- 
mogeneous Bose gas, Tc, with interactions has had a long 
and controversial history jlj-p^ (and further references 
therein). The correction to the ideal gas critical temper- 
ature, T°, is now established theoretically |^,^ (see also 
p| ) to be linear in the scattering length a in leading order 
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where n is the particle density. The variation of Tc is de- 
termined primary by large distance, or small momenta, 
contributions. Due to infrared divergencies around the 
transition point, pertubation theory is not applicable, 
and the coefhcent c cannot be obtained by simple pertur- 
bative techniques. Except in special cases, e.g., a large 
number of internal degrees of freedom, the constant c 
must be evaluated numerically. 

However, even numerical calculations have not settled 
this issue. Whereas Ref. provides c = 2.33 ±0.25 from 
a direct calculation of the coefficient first taking the limit 
ani/^ 0, Ref. jj] predicts c = 0.34 ±0.06 after numeri- 
cal extrapolation of the calculation to the limit a — > 0. In 
this paper we show that the difference in these two results 
is attributable to a non-analytic structure, ^ a^lna, of 
the transition temperature in a. This correction, which 
is negative, does not introduce a new length scale beyond 
n~i/^, but because of its logarithmic character it gives 
rise to a strong dependence on a even in the very di- 
lute limit, av}/^ 0. Very recently, Refs. 11 1 obtain 
c ~ 1.3 from a classical 0^ model on a lattice extrapo- 
lated to the continuum, a result qualitatively consistent 
with the larger value found in Ref. In these calcu- 
lations the logarithmic terms enter as corrections to the 
continuum limit. 



Equation (|I|) is only the beginning of an asymptotic 
expansion, as one might suspect, since for a < the sys- 
tem is unstable. Inclusion of the a? In a term, when one 
extrapolates numerical data from finite a values to the 
limit a — *■ 0, provides a first resolution of the apparent 
discrepancies between numerical calculations done at fi- 
nite a iQ and those valid for a — > ||Jl^,|ll|. To obtain 
a quantitative estimate, we explicitly calculate the log- 
arithmic correction in a model with N internal degrees 
of freedom, to leading order in 1/N . The result suggests 
that the linear increase of ATc at small but finite a is 
noticeably suppressed for the physical case oi N = 2. 

We consider a uniform system of bosons of mass m at 
temperature T , and assume that the two-body interac- 
tion can be described by the s-wave scattering length a. 
Above the critical temperature, the density n is given in 
terms of a sum over Mastubara frequencies, = 2T:ivT 
{v — ±1,±2, . . .), of the single particle Green's function 
G{k,z): 
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where {Ti — ks — 1) 

G~^{k,z) = z + 11 
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and /X is the chemical potential; the condition ^ = S(0, 0) 
determines the Bose-Einstein condensation point. 

The shift of the critical temperature at fixed density 
is more conveniently calculated in terms of the shift, 
Aric — nc{a,Tc) — nc{0,Tc), in the critical density at 
fixed T; the two shifts are related, to the orders of in- 
terest (less than a^), by ATc/Tc = -{2/3)Anc/nc. As 
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shown in |Q, the leading Hnear shift ATc'^'' is given solely 
by the zero Matsubara frequency term: 
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where A = {2tt /mTY/'^ is the thermal wavelength, 
C(3/2) — 2.612 . . ., Go is the Green's function of the ideal 
Bose gas at r°, and 



U{k) = 2m(E(fc,0) - S(0,0)). 



(6) 



At the critical temperature, U{k) can be calculated to 
order a^/A^ by considering only the v = sector, which 
corresponds to a classical field theory. At the transition, 
U{k) has the scaling structure. 



U{k) = ^aikXya), 
(1) 
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from which the linearity of ATc in a follows 
If U{k) is calculated by classical field theory, ATc is 
strictly linear in a. 

(2) 

The next-to- leading order corrections, ATc , arise in 
terms with non-zero Matsubara frequencies, both explic- 
itly [Eq. (j|)] and in internal loops, in the calculation of 
U{k). As we show below, the internal loop corrections 
begin at order and a^lna; the a^lna terms may be 
extracted from 
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[G(fc,z,)-Go(fc,z.)]. (8) 



Since the infrared behavior is regular for v ^ 0, we ex- 
pand the denominator of G to first order in S(fc, z^) — /x, 
and write 



AT 



(2) 



T, 



d^k S(fc,z^)-E(fc,0) + [/(fc)/2m 



iz^-k^/2my 



(9) 



From pertubation theory we know the functional form of 
U{k) outside the critical region in k, 
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this logarithmic behavior is valid to all orders of per- 
tubation theory, as one can verify by power counting. 
Since the dominant contribution to the integral in Eq. (|^) 
comes from momenta k ~ A^^, the ultraviolet behavior 
of U{k) generates a logarithmic shift in the critical tem- 
perature: 
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Since z^) — 0) tends to zero for large momenta 
k, the contribution of this term in Eq. remains of or- 
der a} j}? . Thus the next-to-leading order to the critical 
temperature shift is proportional to a} In a and is always 
negative for small a. 

In order to estimate the shift quantitatively, we calcu- 
late it in the large TV model. The j/ = sector, equivalent 
to a classical (\)^ field theory in three spatial dimensions, 
is described by the action |0| , 
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with u = 96TT^a/X^. The classical field theory suffers 
from ultraviolet divergencies, which can be regularized 
by introducing a large momentum cutoff A. As shown 
in Refs. |@,||, the leading order corrections to the crit- 
ical density are dominated by long distance properties 
and U{k) is independent of the cutoff. Therefore one 
can derive U (k) with a fixed cutoff A in the action, take 
the limit A ^ 00, and determine the corrections to the 
critical density from Eqs. (^) and (||). 

Instead of this procedure we will obtain the next-to- 
leading order corrections in an independent way, which 
has the advantage of making contact with the numeri- 
cal 0** lattice calculations. Starting from the finite tem- 
perature quantum field action, one can derive the ef- 
fective action of the classical field theory by integrat- 
ing pertubatively over the non-zero frequency quantum 
modes, v 0, which provides a large momentum cutoff 
A ^ \J mT ~ 1/A and renormalized effective coefficients 
of the Euclidean action ||l^ . Following |^ , the corrections 
to the transition temperature are given in this effective 
field theory by 
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where the subscript A indicates the explicit dependence 
on the ultraviolet cutoff which incorporates the leading 
effects of non-zero Matsubara frequencies. 

In the large N limit U\{k) is given in terms of the 
particle-hole bubble, B{q), by 
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to leading order in 1/N in three dimensions ul 



1 



6 



8(7 Ng*K 



+ 0(A-2) 



(15) 



where g* = A&tt'^ /N. Following Ref. 
critical temperature shift 
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As we see from the A-dependence of the bubble, Eq. (|15D , 
the A-dependence of U^ik) gives rise to higher order cor- 
rections in A~^. These terms arise effectively from the 
internal non-zero Matsubara frequencies, and lead to cor- 
rections ^ A^^lnA in /^Tc/Tc, which we can neglect. 
Thus 
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where r = (A^u/48)^^. To obtain the leading order cor- 
rections we differentiate with respect to At, which allows 
us to isolate the contributions around Ar oo. Integrat- 
ing back we find 
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The integration constant is given by the large N result 
for the linear shift, Ref. ||^. In contrast to the linear cor- 
rections in a, the next-order terms depend on the short 
length scale properties of the system, modelled by the 
cutoff A, and therefore are not universal. 

To estimate the influence of the logarithmic terms we 
take TV = 2 and AA = (27r)i/2 in Eq. (|l8|). The cor- 
rections to the linear behavior of T, found in this way 
are precisely those one finds for TV = 2 by including the 
particle-hole bubble sum in U{k) . The resulting depen- 
dence of the transition temperature on a is shown in Fig. 
1. For a gas parameter na^ 10~^, corresponding to 
the experimental region of Bose-Einstein condensation 
in atomic gases and liquid ^He in Vycor Jist , and 
to the lowest density Monte Carlo data of Ref. l4|, the 
nonlinear corrections depress the linear shift by ^ 50%; 
instead of the a ^ result c ~ 2.33 in Eq. (|^), one 
obtains a coefficient c ~ 1.2 for an^^^ ~ 10~^. Even if 
the extrapolation from the large N expansion to = 2 is 
unjustified, this calculation suggests that the logarithmic 
terms play an important role for present numerical and 
experimental parameters. The noticeable depression of 
AT, in this parameter regime is also confirmed by self- 
consistent numerical model calculations [O. 



A qualitatively similar strong dependence on a is found 
in the renormalization group calculations of Ref. 
which derives an a In a correction jl^]. Such a result 
would follow from Eq. (||) were U{k) to be linear in 
k up to an ultraviolet cutoff'. However, in the regime 
a/A^ <C A, corresponding to the dilute limit, a/ A <C 1, 
U{k) is given by the perturbative result ( px| ) for momenta 
in the region a/X^ ^ A: ^ A rather than being linear in 
k. 
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FIG. 1. Dependence of the transition temperature of a di- 
lute homogeneous Bose gas on scattering length, from Eq. jig ) 
for N ^ 2. The falloff of ATc/an^^^ with increasing anF^ 
arises from the non-universal next-to-leading order logarith- 
mic corrections. In the regime where the curve is represented 
by a dashed line, higher order corrections begin to become 
important. The dark circle is the calculation of Ref. [6]; the 
points shown as open circles are the lowest density data from 
the numerical calculations of Ref. [4]. The data point shown 
by the cross indicates the numerical results of Refs. [10,11] for 
a lattice theory extrapolated to the continuum. 

The next-to-leading order corrections are important for 
classical 0'* calculations as well p^ , [rT|| ; in the universal 
region where the influence of the ultraviolet cutoff A is 
unimportant (A — > oo), classical field theory provides 
perfect scaling, implying a linear shift of T, for all a. 
Variations in a of AT, in ]To| | arc due to non-universal 
corrections and are sensitive to the details of the scheme 
used to regularize the classical cj)^ theory . Extrapola- 
tion to the universal small coupling region, a/A^A —> 0, 
together with finite size scaling to the thermodynamic 
limit allows one not only to extract the coefficient c of the 
linear shift in T,, but should also provide the magnitude 
of the non-universal a^ In a corrections for the physical 
case N = 2. 
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